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Abstract 

In the framework of the Constrained Minimal Supersymmetric Standard Model 
(CMSSM) we discuss the impact of the recent experimental information, especially 
from E821 Brookhaven experiment on g^ — 2 along with the light Higgs boson mass 
bound from LEP, to the Dark Matter direct searches. Imposing these experimen- 
tal bounds, the maximum value of the spin-independent neutralino-nucleon cross 
section turns out to be of the order of 10~® pb for large values of tan (3 and low 
Mi/2,mQ. The effect of the recent experimental bounds is to decrease the maxi- 
mum value of the cross section by about an order of magnitude, demanding the 
analogous sensitivity from the direct Dark Matter detection experiments. 



Supersyminetry, or fermion-boson symmetry, is an omnipotent and ubiquitous ele- 
ment in our efforts to construct a unified theory of all fundamental interactions observed 
in nature. At very high energies, close to the Planck scale (Mp) it is indispensable in con- 
structing consistent string theories, thus dubbed superstrings. At low energies (~ 1 TeV) 
it seems unavoidable if the gauge hierarchy problem is to be resolved. Such a resolution 
provides a measure of the supersymmetry breaking scale Msusy ~ 0(1 TeV). There is, 
albeit circumstantial or indirect, evidence for such a low-energy supersymmetry breaking 
scale, from the unification of the gauge couplings ^ and from the apparent lightness 
of the Higgs boson as determined from precise electroweak measurements, mainly at 
LEP 0. Furthermore, such a low energy SUSY breaking scale is also favored cosmo- 
logically. As is well known, i?-parity conserving SUSY models, contain in the sparticle 
spectrum a stable, neutral particle, identifiable with the lightest neutralino (x), referred 
as the LSP [0]. One can then readily show ^ that such a LSP with mass, as low-energy 
SUSY entails, in the 100 GeV — 1 TeV region, may indeed provide the right form and 
amount of the highly desirable astrophysically and cosmologically Dark Matter (DM). 
As times goes by, the experimental evidence for DM, from different quarters, strengthens 
in such a way, that it has assumed a central role in the modern cosmology. The most 
recent evidence, coming from the observation of the first three acoustic peaks in the Cos- 
mic Microwave Background (CMB) radiation small angle {6 < (9(1°) anisotropics 0], 
is of tantalizing importance. It is not only provides strong support to a flat [k = 
or fio = 1)5 inflationary Universe, but it also gives an unprecedented determination of 
hdhQ 



^Aihn ~ 0.15 ± 0.05, which taking into account the simultaneously determined baryon 



density fi^/iQ ~ 0.02, and the rather minute neutrino density suggests 

nDMhl = 0.13 ±0.05 (1) 

One then is tempted to combine this recently determined DM density, assuming, as 
we do here, that it is all due to neutralini (i.e. ^dm = ^x)' with other presently 
available constraints from particle physics, in order to find out what is the chances of 
observing, soon or in the near future, DM directly in the laboratory by elastic neutralino- 
nucleus scattering, from the energy deposition in the detectors [Q . These particle physics 
constraints include the lower bound on the mass of the Higgs bosons {rrih > 113.5 GeV) 
provided by LEP [0], the allowed region for b -^ 57, at 95% CL range (2.33 x 10^^ < 



B{b -^ 57) < 4.45 X 10"'^) [0, and the recent results from the BNL E821 experiment [Q 
on the anomalous magnetic moment of the muon (5a^ = 43(16) x 10~^°), assuming, as 
we do here, that is all due (at the 1 or 2 a level) to low-energy supersymmetry. It should 
be stressed that the possibility of a rather sizeable positive contribution to g^ — 2 from 
low energy SUSY in the region of large tan/? and fi > where the b —^ s'y constraint 
weakens considerably, has been long strongly emphasized [0. It is amusing to notice, 
that in our previous analysis of the direct DM searches |jTO|, done before the the BNL 
E821 announcement, we had paid particular interest in the large tan/3 region, since it 
provided the higher possible rates for direct DM detection! Similar results are presented 



in Ref. |]TT|. Actually, as we have stressed for some time now |T^,|T2|, one way to get the 



"right" amount of the neutralino mass density {Q^, h^), even for relative high values of rriQ 
and Ml/2, is to move to the large tan/5 region, because efficient neutralino annihilation 
directly through A and H poles, occurs. The annihilation cross sections increase with 
tan/?: couplings Axx and Arf, Abb increase, while ttla decreases, thus one may also 
expect a rather appreciable increase in the elastic x-^ucleon cross section, as is indeed 
the case. It may turn out, if the BNL E821 result is due to low energy SUSY, that 
the imposed lower bounds on rriQ, M1/2 and lower bounds on tan/? |T^-|T^ make the 



direct neutralino annihilation, through the A, H poles, the major mechanism for getting 



the right amount of DM ||T0|, 0, |T7| , as well as as being consistent with all available 
constraints ||13|,[T^ . 

Before presenting our results we shall give a brief account on the numerical analysis 
employed in this paper. This will be useful in comparing our results with those of other 
authors. In our analysis we use two- loop renormalization group equations (RGE), in the 



DR scheme, for all masses and couplings involved, defining the unification scale Mqut 
as the point at which the gauge couplings di and 02 meet. We do not enforce unification 



of ds gauge coupling with with di^2 at Mqut- The experimental value of the MS strong 
coupling constant at M^, which we consider as input, is related to a^ through as{Mz) = 
a3{Mz)/{l — Ads). Ads represent the threshold corrections which affect significantly 
the value of ds at Mz and hence, through RGE, its value at Mqut- The latter turns 
out to be different from di,2 (Mgc/t), reflecting the fact that gauge coupling uniflcation is 
impossible to implement in the constrained scenario with universal boundary conditions 



for the soft masses. For the determination of the gauge couphngs di,2 we use as inputs 
the electromagnetic couphng constant ao the value of the Fermi coupling constant Gp-, 
and the Z-boson mass Mz- From these we determine the weak mixing angle, through 
s^c^ = Tiao / ^/2M'^G pi^ — Af), and the value of the electromagnetic coupling constant 



at Mz- With Af we denote the correction to the muon decay amplitude. In the DR 
scheme the latter is related to Qq through d(M^) = ao/(l — Adgm), where Adgm are the 
appropriate threshold corrections (see Ref. [|l8l). The input value of the strong coupling 
constant is taken within the experimental range as{Mz) = 0.1185 ± 0.002. 

In running the RGB's, as arbitrary parameters we take are as usual the soft SUSY 
breaking parameters mo, M1/2, ^0 the value of tan/5 and the sign of the Higgsino mixing 



parameter /i. The top and tau physical masses, Mt, M^, as well as the MS bottom run- 
ning mass mb{mh) are also inputs. As default values we consider Mt = 175 GeV, M^ = 
1.777 GeV and mh{mf,) = 4.25 GeV although we allow for variations within their exper- 
imentally allowed region. The determination of the bottom and tau running masses at 



Mz is done by running S'f/c(3) x f/em(l) MS RGB's, using three-loop RGB's for the 
strong coupling constant. We also include two- loop QBD corrections, as well as two-loop 



contributions from the interference of the QCD and QBD corrections. The running MS 



masses at Mz are then converted to DR in the usual way. From these we can extract 
the corresponding Yukawa couplings at Mz- We point out that the important QCD as 
well as the supersymmetric gluino, sbottom and chargino, stop corrections to the bottom 
mass are duly taken into account. For the determination of the top Yukawa coupling 
at Mt we relate its pole and running masses taking into account all dominant radiative 
corrections. By running the RGB's we can have the value of the top Yukawa coupling 
at Mz- 

The determination of the Higgs and Higgsino mixing parameters, rn^ and /i, is a 
more subtle issue. These are obtained by solving the minimization conditions with the 
one- loop corrected effective potential with all particle contributions taken into account. 
Since large values of tan f3 cause large logarithmic corrections, invalidating perturbation 
expansion, we solve the minimization equations taking as reference scale the average 
stop scale Qt — y/fnfjri^. At this scale the corrections are numerically small and hence 
perturbatively valid. Thus in each run we determine m\{Qi), ^{Qi)- The values of 



mllMz), fi{Mz), whenever needed, can be found by solving the RGE's having as initial 
conditions the values of these quantities at Q^. 

For the calculation of the lightest supersymmetric particle (LSP) relic abundance, we 
solve the Boltzmann equation numerically using the machinery outlined in Ref. [|T2|. In 



this calculation the coannihilation effects, in regions where r^ approaches in mass the 
LSP, which is a high purity Bino, are properly taken into account. 

Before embarking to analyse our numerical findings it would be beneficial to review 
the physical mechanism through which the scalar, i.e. spin-independent, x-nucleon cross 
section (agcaiar) is enhanced, to levels approaching the sensitivity of ongoing experiments. 



The (J scalar IS enhanced in the region of the parameter space where tan/3 is large [|I0 
Q. The dominant contribution to this regime is the Higgs boson exchange. For given 
inputs mo, Ml/2, ^0 and the sign of /i, Higgs masses decrease as tan/5 increases. Hence 
the contribution of Higgs bosons to neutralino-quark elastic cross section becomes more 
important in the large tan (3 regime. Such a decrease in the mass is not sufficient by itself 
to increase (Jgcaiar- The major role in this increase plays the coupling of the CP-even 
heavy Higgs whose coupling to rf-quark is proportional to ^^, which is proportional to 
tan/3, when the latter becomes large. The coupling of the light CP-even Higgs, unlike 
the heavy Higgs case, does not grow with increasing tan/3 but stays constant of order 
unity. Therefore despite the fact that the heavy CP-even Higgs is heavier than its light 
CP-even counterpart, its contribution may be much larger in the large tan/3 region, due 
to its enhanced coupling to rf-quark. In the large tan/3 regime the mass of the heavy 
CP-even Higgs can be approximated W^ by the relation 



niH^m^ + mA-mz-e, (2) 

where e represents the leading stop 1-loop corrections to the CP-even Higgs masses. From 
this is obvious that the lowest tuh values are obtained in the region where itia is light and 
rrih is close to its lower experimental bound. Moreover it is worth pointing out that in the 
large tan/3 region neutralino relic densities decrease as we have already emphasized due 
to both the decrease of the pseudoscalar mass, whose exchange in xx ^ bb , rf processes 



^Enhancement of a scalar is also possible in the context of the so-called focus point supersymnietry 
scenario iQ, where toq > 1.5 TeV, yet such large values of ttiq are not favourable by the recent g^ — 2 
data. 



is less suppressive, and the increase of the XX^ as well as the Abb and Arf couplings. 
The smallness of the LSP's Higgsino component is compensated by the largeness of tan/3 
yielding neutralino annihilation cross sections compatible with the recent astrophysical 
data. Hence there are regions in which we can obtain both low relic densities and high 

^scalar • 

Bearing all these in mind, we proceed discussing our findings. For our numerical 
analysis a large random sample of 45,000 points in the region of the parameter space 
designed by 2 < tan/3 < 55, M1/2 < 1.5 TeV, mo < 1.5 TeV, \Ao\ < 1 TeV, and /i > 
is used. The /i < case is not favored by the recent b —* s'j data, as well as by the 
observed discrepancy of the g^ — 2, if the latter is attributed to supersymmetry, and 
therefore we shall not discuss it in the sequel. It is also worth noticing that in the /i > 
case the constraint from 6 — > 37 data is superseded by the ruh > 113.5 GeV bound, 
in the bulk of the parameter space |13[ . In figure [I| we plot the scalar x-nucleon cross 



section as function of the LSP mass, rriy.. On the top of the figure the shaded region (in 
cyan colour) is excluded by the CDMS experiment |^. The DAMA sensitivity region 
(coloured in yellow) is also plotted [^. Pluses (+) (in blue colour) represent points 
which are both compatible with the E821 data aj^^^ = (43.0 ± 16.0) x 10"^° and the 
cosmological bounds for the neutralino relic density ^x^o ~ 0.13 ± 0.05. Diamonds 
(o) (in green colour) are points which are cosmologically acceptable with respect to the 
aforesaid bounds, but the bound to the a^^^"^ has been relaxed to its 2a region, namely 
11 < a^^^^ X 10^° < 75. The crosses (x) (in red colour) represent the rest of the points 
of our random sample. Here the Higgs boson mass bound, rrih > 113.5 GeV has been 
properly taken into account. From this figure it is seen that the the points which are 
compatible both the g^ — 2 E821 and the cosmological data (crosses) yield cross sections 
of the order of 10^^ — 10^^ pb and the maximum value of the m^ is about 200 GeV. If 
one considers the 2a region of the g^j, — 2 bound the preferred cross sections can be as 
small as 10^^° pb and correspondingly the upper bound of m^ is drifted up to 350 GeV. 
In the following figures |^ and ^ the ascaiar is plotted as function of the parameters mo 
and tan/3 respectively. One can see that the points which conform to cosmological and 
la g^ — 2 experimental constraints, yield a maximum value of itiq about 600 GeV, and 
for the 2a case 1200 GeV. The aforementioned bounds on the m^ and uiq are related 



with the analogous bounds put on the soft parameter Mi 12 and mo from the g^ — 2 E821 



data pi 0,116 



From figure ^ it is apparent that the majority of the points that are compatible 
with cosmological and g^ — 2 data are accumulated toward rather large values of tan /5, 
specifically tan/? > 40, although there are indeed few of them with smaller values of tan /3. 
As it has been already pointed out in the large tan (3 region we can have simultaneously 
cosmologically acceptable values of fi-^ h^ and also big values for the elastic cross section 
X-nucleon. Furthermore the g^ — 2 muon data prefer large values of tan/5, as «^^^^ 
is proportional to tan/? [^]. Therefore as tan/? increases large regions of the parameter 
space (mo, M1/2) are compatible with the E821 experimental constraints. Taking all these 
into account it is not surprising that the conjunction of the cosmological and g^ — 2 data 
happens for large values of the tan (3 and for large scalar cross section x-nucleon, as it 
can be perceived from figure |^. 

Comparing figure |I| and ^ one can realise how g^ — 2 data constrain m^^ mass to be 
up to 200 GeV or 350 GeV for the Icr or 2a case respectively. In figure ^ we don't 
impose the constraints stemming from g^ — 2 data, therefore due to the coannihilation 
processes the cosmologically acceptable LSP mass can be heavier than 500 GeV. What 
is also important to be noticed about the direct searches of DM is that imposing the 
g^ — 2 data the lowest allowed x-nucleon cross section increased by about 2 orders of 
magnitude, from 10^^^ pb to 10^^ pb. This fact is very encouraging for the future DM 
direct detection experiments. Figure ^ illustrates the significance of the Higgs boson 
mass bound. If one allows for values m/j > 100 GeV many points which yielding cross 
sections even (9(10~^) pb appear. Comparing figure |I] and ^ we observe that the recent 
Higgs mass bound [mh > 113.5 GeV) reduces the maximum value of the scalar cross 
section for about one order of magnitude, that is from 10"^ pb to 10~^ pb. There is 
direct and indirect dependence of the a scalar on rrih. As nih decreases its contribution to 
the ascaiar, being proportional to 1/m/j, increases leading to larger agcaiar- The indirect 
relation of the ascaiar to mh can be perceived from Eq. ^ Light m^, results to light niH 
and therefore to large ascaiar again. 

Concluding we have studied the impact of the recent experimental information to the 
DM direct searches. Especially we have considered the effect of the recently reported 



deviation oi g^ — 2 from its SM value, as well as of the light Higgs boson mass bound 
from LEP experiments. The imposition of these experimental constraints results to a 
maximum value for the spin-independent X'^iucleon cross section of the order of 10^^ 
pb for uiy. ~ 100 GeV as small as allowed by chargino searches, and for tan/3 > 45 as 
large as possible for the Higgs states to be as light as allowed by theoretical constraints 
and experimental searches. As it can be seen from figures § and |^ the effect of these 
experimental constraints is to decrease the maximum value of (Tscaiar by almost one 
order of magnitude, and therefore to make the direct detection of the LSP on the future 
experiments by some means more difficult. 
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Figure 1: Scatter plot of the scalar neutralino-nucleon cross section versus rriy., from a 
random sample of 45,000 points. On the top of the figure the CDMS excluded region 
and the DAMA sensitivity region are illustrated. Pluses (+) are points within the E821 
experimental region a^^^"^ = (43.0 ± 16.0) x 10^^'^ and also cosmologically acceptable 
r2^ /iq = 0.13 ± 0.05. Diamonds (o) are also cosmologically acceptable points, but with 
Q,susY ^i^j^in lY^Q region 11 x 10~^° < a^^^^ < 75 x 10"^°. Crosses (x) represent the 
rest of the random sample. The Higgs boson mass bound ruh > 113.5 GeV is properly 
taking into account. 
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Figure 2: In this figure we display the scalar neutralino-nucleon cross section versus mo- 
Points are as in Fig. |l|. 
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Figure 3: In this figure we display the scalar neutralino-nucleon cross section versus 
tan/?. Points are as in Fig. |l|. 
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Figure 4: Scatter plot of the scalar neutralino-nucleon cross section versus m^, from a 
random sample of Fig. |I]. Diamonds (o) are cosmologically acceptable points, without 
putting an restriction from the a^^^"^. Crosses (x) represent points with unacceptable 
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Figure 5: Scalar neutralino-nucleon cross section versus m^. Points are as in Fig. |T[ 
Here the Higgs boson mass bound {rrih > 113.5 GeV) has been relaxed and the bound 
rrih > 100 GeV is used. 
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